arXiv:1506.03747vl [physics.atom-ph] llJun2015 


Multi-electron dynamics in the tunnel ionization of correlated 

quantum systems 

Maximilian Hollstein and Daniela Pfannkuche 
Universitdt Hamburg, I. Institut fiir Theoretische Physik, 

Jungiusstrafie 9 , 20355 Hamburg, Germany 
(Dated: June 12, 2015) 

Abstract 

The importance of multi-electron dynamics during the tunnel ionization of a correlated quantum 
system is investigated. By comparison of the solution of the time-dependent Schrodinger equation 
(TDSE) with the time-dependent configuration interaction singles approach (TDCIS), we demon¬ 
strate the importance of a multi-electron description of the tunnel ionization process especially for 
weakly confined quantum systems. Within this context, we observe that adiabatic driving by an 
intense light field can even enhance the correlations between still trapped electrons. 
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INTRODUCTION 


The ionization of closed shell atoms is impressively well-understood on the basis of sin¬ 
gle active electron approaches (SAElJ-yl) or effective one-particle theories as the time- 
dependent conhguration interaction singles approach (TDCIsQ-Q). Within these ap¬ 
proaches, ionization is described by the ejection of a single electron into the continuum 
while the residual electrons remain unaffected (i.e. they are only taken into account by 
a time-independent potential for the active electron (SAE) or they are kept residing in 
Hartree-Fock ground state orbitals of the held-free atom (TDCIS)). However, in weakly 
conhned quantum systems such as molecules or atom-like systems as semiconductor quan¬ 


tum dots in which the e 


the trapped electrons 


ectron-electron interaction induces signihcant correlations between 


llj, an independent particle description of the ionization process. 


as inherent in these approaches, is expectably unsuitable from the very beginning. In this 
paper, we address the question concerning the importance of a multi-electron description 
of the tunnel ionization process of weakly conhned and correlated quantum systems. For 
this purpose, we consider the dynamics which is induced by an intense low-frequency light 
held in a two-electron model system. Our conclusions drawn from these considerations, 
however, are not only valid for two-electron quantum systems but they are transferable to 
systems with more than two weakly bound electrons. Furthermore, we would like to remark 
that the conclusions drawn from the presented model calculations are not limited to the 
ionization of molecular or atomic systems but they are also relevant for experiments basec 


on the application of a ('time-dependent) yoltage to semiconductor quantum dots 


and photo assisted tunneling (PAT) 
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, Il4l-ll7l| in the low-frequency regime. 


This paper is structured as follows: after introducing the considered model system, we 
demonstrate the need for a multi-electron description of the tunnel ionization process by 
comparison of the solution of the time-dependent Schrodinger equation (TDSE) with the 
results obtained by the TDCIS approach. By detailed analysis of the exact wavefunction we 
provide insight into the light induced dynamics revealing that a multi-electron description 
becomes necessary not only due to ground state correlations but also due to a light in¬ 
duced collective electron motion which is accompanied with an enhancement of correlations 
between still trapped electrons. 
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MODEL 


In order to study the ionization of a weakly confined quantum system, we consider a 
one-dimensional model system consisting of two electrons in an inverse Gaussian confining 
potential. (Note that the effective potential realized in semiconductor quantum dots can 


be well approximated by this potential [l8|, ll9|). Without exception, atomic units are used 


throughout this paper. Referring to the spatial coordinate as x, the considered conhning 
potential can be denoted by: 


V{x) = -Voe 




( 1 ) 


where Vq determines the depth and w the width of the conhnement. Within the harmonic 
approximation of this potential, the strength of the conhnement i.e. the level spacing of the 
lowest bound states is given by: 

AE ^ ( 2 ) 


w 

In the subsequent considerations we make use of this dependency and vary the strength of 
the conhnement by variation of the width w. 

The electron-electron interaction is taken into account by a regularized Coulomb interaction: 

1 


Vee{Xl,X2) = 


( 3 ) 


\/{xi - X2y + 5^ 

where xi ^2 denote the spatial coordinates of the two electrons. By choosing 6 unequal zero 
(we set 6 to 0.5), a hnite width of the electronic wave function along the unconsidered 
spatial dimensions is taken into account in a phenomenological wavj^. 211. For the results 
presented below, the depth of the conhning potential is set to Vq = 3. 


Hereinafter, we restrict our considerations to light helds with a wavelength which is large 
compared to the extend of the considered quantum system. This allows the application 
of the dipole approximation to the light-matter interaction. Thus, in length gauge, the 
Hamiltonian of the quantum system interacting with the light held is given by: 




2^dx^ 

2=1 '' 




( 4 ) 


( 5 ) 


2=1 
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where Hq is the Hamiltonian of the held free system and E(t) denotes the time-dependent 
electric held of the laser. Subsequently, we consider the initial dynamics of the system in an 
electric held with sinusoidal time dependence i.e.: 

E{t) = i?osin(a;f). (6) 


In the following, we consider the situation where the two-electron system is initially prepared 
in the singlet ground state. In order to study the light induced dynamics, we solve the time- 
dependent Schrodinger equation numerically. For this, we employ a split operator technique 


22 ( 1 . Rehections at the boundaries of the hnite simulation box are avoided by the use of a 


complex absorbing potential Vcap- 


-iC(||x|-xol)^ \x\>Xq 

Vcap — < (7) 

I 0 else 

For the calculations presented below, C is chosen to 4 x 10“^. Xq is chosen to 15 where both 
the conhning potential and the ground state wavefunction are only insignihcant diherent 
from 0 (i.e. for all model parameters, the absolute value of conhning potential is for x = 15 
smaller than 0.001 and the value of the ground state density is decreased to ^ 10“^® ). 

In this paper, we consider the tunnel ionization process i.e. the ionization of the considered 
quantum system via tunneling through the ehective potential barrier which is formed by 
the instantaneous electric held and the conhning potential. Considering that the tunnel 
ionization regime is characterized by a Keldysh parameter much smaller than one, we 
choose the frequency of the light held to ca = 27r x 0.001 and the amplitude to Eq = 0.3 
corresponding to a Keldysh parameter smaller than 0.01 for the considered model potential 
parameters. For semiconductor quantum dots, the characteristic conhnement energy is in 
the meV range so that the considered light held would be in the domain of far infrared light to 
microwaves. We found that then, excitations within the potential well are negligible for the 
considered laser parameters so that we here consider subsequently only the quasi-adiabatic 
dynamics during a half-cycle of the laser held. Since for tunnel ionization, the ionization 
process is most prominent at times when the electric held of the light held is extremal, we 
focus especially on the situation at t = T/4 when the electric held is maximum. Here, we 
restrict our considerations to the dynamics in the conhning potential well and its vicinity 
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by analyzing the restricted wavefunction: 




d'(Xi,X2,t)||xi,2|<x„ 


/jx j. dx,dX2j^(Xr,X2)j^ 

y ^ ‘^max ^ •^max 


whereas nnless noted differently, we chose 


^max 


= Xq = 15. 


( 8 ) 


( 9 ) 


Before the two-electron system is exposed to the light field, both electrons are sitnated in 
this area. However, at later times, electron(s) are excited by the laser into the continuum 
and leave the considered spatial region. We regard these electrons as ionized and determine 
the ionization probability consequently by: 

/ ^max r^max 

/ dxidx2\^{xi,X2,t)\^. ( 10 ) 

'^max ^max 


RESULTS 

Comparison of TDCIS and TDSE 

Within the TDCIS approach, the two-electron wavefunction is expanded in the Hartree- 
Fock ground state |<ho) and its particle-hole excitations |<Fo) 

|>I>>=a„H>„>+5^aiH>-> (11) 

i>0 

with 

l^o) = ^(4tCot + clwoi) l^o) (12) 

where cof, cqj, denote annihilation operators of the spin orbitals which are occupied in 
the Hartree-Fock ground state determinant while c|.j., denote creation operators of virtual 
orbitals. The spatial part of the TDCIS singlet two-electron wavefunction is consequently 
given by: 

T(a;i,a;2,t)j,oc/s = ['0o(u)'0i(a^2) + (13) 
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with 


ai{t) = 


(14) 


ao{t)/2 i = 0 
a^it) else. 

As can be seen in equation ffT^ . at all times, only one electron can be active (i.e. occupy 
an arbitrary orbital) while the other electron is forced to occupy the Hartree-Fock ground 
state orbital -^o- 

In order to determine the coefficients ai{t) we solve the Hartee-Fock equations on a pseudo- 
spectral grid as described in [2^. The TDCIS wavefunction is then propagated by iterative 


Lanczos reduction 


25| within the configuration interaction singles singlet subspace con¬ 


structed by the eigenfunctions of the Fock-operator. 

In hgure [H the ionization probability as defined in eq. [lO] obtained by the exact solution of 
the TDSE and by the TDCIS approach are shown for varying confinement widths. As can 
be seen, the TDCIS approach reproduces very accurately the exact ionization probability for 
w < 2.5 whereas large deviations are observable for w > 2.5 (e.g. by 40 percent for w = 5). 



FIG. 1. (Color online) The ionization probability Pion as defined in equation 1101 after a half cycle 
of the laser field obtained by the solution of the TDSE (red curve) and by the TDCIS approach 
(blue curve with dots) in dependence on the width of the confinement. For narrow confinements 
(i.e. for w < 2.5), the TDCIS approach reproduces very accurately the exact ionization probability 
whereas for wide confinements substantial deviations are observable. 

Some insight can be gained by the approximation of the exact two-electron wavefunction 
in terms of conhgurations constructed from the two most occupied natural orbitals (cf. eq. 
fTS]) . This truncated conhguration interaction expansion allows a good approximation of the 
exact wavefunction since the naturals orbitals constitute an orbital basis set leading to the 
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most rapidly converging expansion in configurations 261 


\1/ ~ E Cij(t)i{xi)(t)j{x2)/ 

ij=l 


\ 


E 


I 


( 15 ) 


If the coefficients c,j are real, this approximate wavefunction can be represented exactly 
in open shell form 27j where one electron occupies a spatial orbital while the other 
electron occupies a different orbital <py, which is not nessecary completely orthogonal or 
totally parallel to 0^: 


'h(xi,X2,t) ^ + 4)y{xi,t)4)u{x2,t). 


(16) 


Although the coefficients Cij of the approximate wavefunction (lT5l) obtained from the nu¬ 
merical propagation are complex at t = T/4, we found for this situation that with good 
accuracy Cm and Cn can be chosen real with alternating sign. Since Cqi and Cio vanish ex¬ 


actly 


28 


29|). a representation in open shell form is nonetheless possible. As described in 


reference |27t], and can be chosen as: 


0U — (|cooh0o + |cilh0i)/(4(|coo|^ + |ciip))4 (17) 

4>v = (|cooh0o “ |cii| 2 0 i)/(4(|coo|^ + |cii|^))4 . (18) 

With this, both wavefunctions (eq. Upland eq. [TB]) have an overlap of 0.967 with the exact 
wavefunction and thus provide a reasonable approximation. Now, the approximation of 
the exact wavefunction in open-shell form allows a convenient comparison with the TDCIS 
wavefunction which can be exactly represented in open-shell form: 


^{Xi,X2,t)^^^jg = 'll:u{xi,t)'lljy{x2,t) + 'll:y{xi,t)'lljy{x2,t) (19) 

where: 

= "00 (20) 

and 

iJu{t) = ( 21 ) 

i 

The two orbitals <pu and (j)y and respectively and ipy are only hxed up to a factor i.e. a two- 
electron wavefunction in open shell form is invariant under the transformation 0^ 4>u/ci, 

(f)y a(j)y where a is an arbitrary non-zero complex number. Therefore, we compare the 
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FIG. 2. (Color online) The renormalized orbitals which constitute the (approximate) open-shell 
form of the exact wavefunction and the TDCIS wavefunction at t = r/4 for a confining potential 
width w = 5. Both the exact and the TDCIS wavefunction are described by one orbital which is 
localized in the potential well (l)v and respectively "0^ (solid curves) and one orbital which has a 
significant overlap with the tunneling barrier describing an tunneling electron and respectively 
i/ju (dashed curves). Both orbitals obtained from the exact wavefunction are shifted to the right 
with respect to the corresponding orbitals obtained from the TDCIS wavefunction. 

TDCIS wavefunction with the exact wavefunction by consideration of the corresponding 
renormalized orbitals (cf. hgure[2]). As can be observed in hgure [2] for w = 5, within this 
open shell approximation, both the exact and the TDCIS wavefunction are described by one 
orbital which is localized in the potential well (0„(TDSE) and respectively (TDCIS)) and 
one orbital which has a signihcant overlap with the tunneling barrier describing an tunnel¬ 
ing electron (0„(TDSE) and respectively 'ifju (TDCIS)). This indicates that the dominant 
ionization process is also within the exact treatment a process where one electron resides in 
the potential well while the other is tunneling through the potential barrier. For TDCIS, the 
localized orbital (i.e. 'ip^) coincides with the Hartree-Fock ground state orbital 'ipo (cf. eq. 
[20]) which is centered at the origin {{x)^ = = 0) while the corresponding orbital of 

the open shell approximation of the exact wavefunction ((^)<^„ = 0.862) is shifted towards 
the local minimum of the instantaneous potential at a; = 1.344. Signihcant deviations are 
also observable between the orbitals describing the tunneling electron i.e. for the TDCIS 
approach and (j)u for the open shell approximation of the exact wavefunction. In particular 
the orbital 0^ (TDSE) ((x)^ = 4.829) is shifted towards the potential barrier maximum 
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at x = 6.385 with respect to V'm (TDCIS) = 4.2026) resulting in a larger overlap of 

(j)u (TDSE) with the potential barrier. This indicates that the immobility of the localized 
electron within the TDCIS approach effects via the Coulomb interaction also the tunnel¬ 
ing electron. That is, within the TDCIS approach, the tunneling electron is less effectively 
pushed towards the tunneling barrier. This provides an explanation for the considerable 
differences between the ionization probability obtained by the TDCIS approach and the 
exact treatment. That is, in particular the underestimation of the ionization process by the 
TDCIS approach can be related to this circumstance(cf. figure [T]). 

So far only an open shell approximation of the exact wavefunction is considered which ap¬ 
proximates the full wavefunction fairly well (i.e. the overlap with the exact wavefunction is 
0.967) but obviously there are still some deviations to the exact wavefunction. In the follow¬ 
ing section we therefore provide more extensive insight into the light induced multi-electron 
dynamics from a different perspective by considering the complete wavefunction obtained 
from the exact treatment. 


Field-Induced Multi-Electron Motion 


The low-frequency field considered here induces a a quasi-adiabatic electron motion lead¬ 
ing to a time-dependent shift of the electronic center-of-mass. Within the harmonic approx¬ 
imation of the confining potential and within the dipole approximation for the light held, 
this shift is given by: 


Ax 



( 22 ) 


Hence, the amplitude of the collective two-electron motion increases monotonically with 
increasing width of the potential well (i.e. oc tc^). Noteworthy, this shift of the potential 
minimum results in the motion of both electrons which becomes noticeable in the open shell 
approximation (cf. eq. [T^ by the fact that the orbital (p.^, which describes a localized 
electron, is not centered at the origin but is rather shifted towards the local minimum of 
the instantaneous potential (see hg. [2]). The regime of large deviations of the ionization 
probability obtained by TDCIS and exact treatment coincides with the regime where Ax 
is in the order of the Bohr radius or larger i.e. Ax ^ 1 = w ^ 4.5. This supports the 
explanation for the large deviations between the ionization probability obtained by TDCIS 
and exact treatment as given in the previous section. That is, within the TDCIS approach. 
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the tunneling electron is less effectively pushed towards the tunneling barrier due to the 
immobility of the localized electron. 

However, since the considered conhning potential well is not a pure harmonic potential 


also excitations of the relative motion are possible 


18| . This allows for correlated electron 


dynamics. In order to get insight into this light-induced correlation dynamics, we consider 
electronic correlations which become noticeable by the circumstance that the electronic 
wavefunction can not be represented by a single Slater determinant. Therefore, we determine 

the degree of oo.xelat.ou ofp two-eleettou wavefunction by eu.ploy.ug the uteaeute of 
correlations K as defined in ^ 


K = (23) 

i 

where rii denote the occupation numbers of the natural orbitals i.e. the eigenvalues of the 


first order density matrix 


31| . Thereby, K can be interpreted as the ’’number” of Slater 


determinants which are effectively necessary to represent the wavefunction (cf. 3^). Since 
every fully uncorrelated singlet two-electron wavefunction is a Slater determinant with a 
doubly occupied spatial orbital, the absence of correlations is characterized by a measure of 
correlations K equals one. Thus, electronic correlations manifest themselves by a value of K 
larger than one (Note that n* > 0 and = 1 cf. M] and [^). In order to determine K 

numerically, we obtain the occupation numbers of the natural orbitals by diagonalization of 


the hrst order density matrix 


27( 1 represented on a spatial grid. 


In hgure [3l K is shown in dependence on the width w of the confining potential for two 
situations i.e. for the singlet ground state (blue curve) and for the two-electron wavefunction 
at t = T/4 (red curve). As can be seen in figure [3l for narrow confinements (tc ~ 1) the 
measure of correlations K for the singlet ground state is nearly one indicating an accurate 
description by the Hartree-Fock determinant (for w = 1 one finds indeed that mostly only 
one natural orbital is populated in the ground state with an occupation of 0.992 ). However, 
K increases monotonically with the conhning potential width. This is related to a feature 
well known for the harmonic conhnement potential, namely, that the conhnement energy 
and the Coulomb interaction scale differently with respect to the characteristic conhnement 
length Iq. That is, whereas the conhnement energy scales as A, the Coulomb interaction 
is proportional to Hence, the relative strength of the Coulomb interaction increases 
monotonically with the width of the conhnement potential leading to a correspondingly 
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FIG. 3. (Color online) The degree of correlations K as defined in eq. [23] for the singlet ground 
state (blue curve with dots) and at the time when the electric held of the laser is strongest (t = 
T/4) (red curve). Since the relative strength of the Coulomb interaction increases with increasing 
width of the conhning potential, the amount of correlations present in the ground state increases 
correspondingly (cf. blue curve with dots). While for narrow conhnements (tc ~ 1), the amount of 
correlations remains constant, for broad conhnement potentials, in presence of the light held, the 
amount of correlations is signihcantly increased in comparison to the ground state. 


increasing population of more than one natural orbital (cf. figure |4j). For w = 5, one 
finds that a second natural orbital is significantly populated with an occupation of 0.107 
so that here, a multi-determinant treatment is already necessary to represent the ground 
state accurately (Note that the overlap with the Slater determinant with maximum overlap 
with the exact wavefunction is given by y/nj = 0.9398 where rii is the occupation number 
of the most occupied natural orbital (cf. j27|) ). However, if one considers the two-electron 
wavefunction at f = T/4, one observes a significant enhancement of correlations with respect 
to the ground state correlations. For instance for w = 5, the occupation of the second most 
occupied natural orbital increases from 0.107 to even 0.142 in the presence of the light 
field and furthermore also a third natural orbital becomes noticeably populated with an 
occupation of ~ 0.028 whereas its ground state occupation is only 0.009 (cf. hgure[5|). Note 
that the occupation of a third natural orbital is not included in the open shell approximation 
(cf. eq. mi) and thus represents an additional deviation between TDCIS wavefunction and 
exact wavefunction which comes along those discussed in the previous section. Since the 
degree of correlations at t = T/4 is independent of the laser frequency for frequencies 
smaller than 27r x 0.0025 (cf. hg. [6|), non-adiabatic excitations within the potential well 
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FIG. 4. (Color online) The probability density p of the singlet ground state and the densities of the 
first three most occupied natural orbitals (j)i, cj }2 and (/)3 weighted by their occupation numbers (ni,n 2 
and 713 ) for a relatively wide confining potential (w = 5). Due to a strong Coulomb interaction 
with respect to the confinement energy, more than one natural orbital is significantly occupied in 
the ground state. 



FIG. 5. (Color online) The probability density p of the two electron state at the time when the 
electric field of the laser is strongest (t = T/4) and the densities of the first three most occupied 
natural orbitals (l)i, 4>2 and (^3 weighted by their occupation numbers for a relative wide confining 
potential (re = 5). In comparison to the ground state (cf. figure |4|), the occupation of the second 
and the third most occupied natural orbital is increased. 

are appearently not the reason for the enhancement of K (the considered laser frequency is 
2t[ X 0.001). Since the most occupied natural orbitals are localized within the well of the 
instantaneous potential, this effect can be attributed to a quasi-adiabatic rearrangement of 
the still trapped electrons within this strongly deformed potential well. The enhancement 
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of the degree of correlations indicates here an effective broadening of the potential well by 
the light held which results in an enhanced effective strength of the Coulomb interaction 
between the still trapped electrons. 



FIG. 6. (Color online) The degree of correlations K (cf. eq. [23]l at t = T/4 for w = 5. Non- 
adiabatic excitations within the potential well manifest themselves by a frequency dependence of 
K for T < 400 i.e. for frequencies larger than lir x 0.0025. Thus, for the considered laser frequency 
(27r X 0.001 - cf. red marker) non-adiabatic effects appear to be irrelevant. 
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CONCLUSIONS 


To conclude, we demonstrated here that for an accurate description of the tunnel ion¬ 
ization process, a multi-electron description is the more important the weaker the conhning 
potential. By analysis of the exact wavefunction, we show that a multi-electron description 
is not only necessary due to ground state correlations but also due to a collective and cor¬ 
related multi-electron motion resulting from a strong deformation of the conhning potential 
well by the light held. 
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